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METHODS FOR IMPROVING 
CONTRAST BASED DYNAMIC 
RANGE MANAGEMENT 



ickground of Invention 

This invention relates generally to contrast-based dynamic range management (C- 
DRM), and more particularly to a method, a related memory medium, and an 
apparatus for implementing and performing a contrast-based dynamic range 
management algorithm. 

Digital x-ray imaging is a well-known, non-contact technique for observing, in 
real time, interior aspects of an object. In practice, an x-ray beam is generated and 
targeted on an object of interest. A detecting device is positioned on the other side of 
the object and detects the x-rays transmitted through the object. The detected x-ray 
signals are converted to digital signals that represent various features in the object, 
are further processed, and the resulting signal is displayed on an image display device 
such as a CRT. 

One of the fundamental image processing problems in digital x-ray imaging is the 
need to transform the intensity dynamic range of the input image to the dynamic 
range supported by an available display device. Typically, the intensity dynamic range 
of an image exceeds the dynamic range of the display mechanism by several times. 
The objective of the transformation accordingly is to compress the "DC" or mean 
component of the different regions comprising the image so that the dynamic range 
(typically 256 gray levels) of an available display device may be utilized in a fashion 
which maximizes the displayed image contrast taking into account, as well, the 
generally recognized limitation of the human eye of discerning only 256 gray levels. 
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[0004] Common approaches to achieving such transformations are known as "unsharp 
masking" and "extended dynamic range (EDR)", the latter a special adaptation of the 
former. The conventional approach to performing the EDR algorithm, in general, is 
simply to subtract a portion of the mean from the input signal. However, if the input 
image is square root encoded, a contrast gain can be added to compensate to some 
extent for the non-linear attenuation of the contrast signal and subsequent contrast 
deficit. 

Summary of Invention 

[0005] In one aspect of the invention a method is provided for compressing an intensity 
dynamic range of an input image to a reduced intensity dynamic range of an image 
display device. The method comprises defining a plurality of units of the input image, 
determining a local mean estimate of an intensity of each of the plurality of units of 
the input image, generating a contrast modification processing value for each local 
mean estimate, and generating an output intensity value, for each local mean 
estimate, by using only the generated contrast modification processing value and an 
input intensity value. 

[0006] In another aspect, an apparatus is provided for compressing an intensity dynamic 
range of an input image to a reduced intensity dynamic range of an image display 
device. The apparatus comprises an image detector, and a computer coupled to the 
image detector. The computer is configured to define a plurality of units of the input 
image, determine a local mean estimate of an intensity of each of the plurality of units 
of the input image, generate a contrast modification processing value for each local 
mean estimate, and generate an output intensity value, for each local mean estimate, 
by using only the generated contrast modification processing value and an input 
intensity value. 

[0007] 

In yet another aspect, a computer readable medium encoded with a program 
executable by a computer for compressing an intensity dynamic range of an input 
image to a reduced intensity dynamic range of an image display device is provided. 
The program is configured to instruct the computer to define a plurality of units of the 
input image, determine a local mean estimate of an intensity of each of the plurality of 
units of the input image, generate a contrast modification processing value for each 
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local mean estimate, and generate an output intensity value, for each local mean 
estimate, by using only the generated contrast modification processing value and an 
input intensity value. 

[0008] In a further aspect, a method for compressing an intensity dynamic range of an 
input image to a reduced intensity dynamic range of an image display device is 
provided. The method comprises defining a plurality of units of the input image, 
determining a local mean estimate of an intensity of each of the plurality of units of 
the input image using a morphological filter, generating a contrast modification 
processing value for each local mean estimate, and generating an output intensity 
value, for each local mean estimate, by using the generated contrast modification 
processing value and an input intensity value. 

[0009] In a still further aspect, a method for compressing an intensity dynamic range of 
an input image to a reduced intensity dynamic range of an image display device is 
provided. The method comprises defining a plurality of units of the input image, 
determining a local mean estimate of an intensity of each of the plurality of units of 
the input image, generating a contrast modification processing value for each local 
mean estimate using a sigmoid look-up table, and generating an output intensity 
value, for each local mean estimate, by using the generated contrast modification 
processing value and an input intensity value. 

Brief Description of Drawings 

[001 0] Figure 1 is schematic flow diagram of an x-ray imaging system implementing and 
performing a C-DRM algorithm in accordance with the invention. 

[001 1 ] Figure 2 is a flow diagram of an EDR algorithm processor of the prior art. 

[001 2] Figure 3 is a flow diagram of a method of implementing and performing the C- 
DRM algorithm in accordance with an exemplary embodiment of the invention. 

[001 3] Figure 4 is a flow diagram of an apparatus implementing and performing the C- 
DRM algorithm, in accordance with an alternative embodiment of the invention. 

[001 4] Figure 5 is a flow diagram of an apparatus implementing and performing the C- 
DRM algorithm, in accordance with an exemplary embodiment of the invention. 
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Detailed Description 



[001 5] Figure 1 is a schematic flow diagram of an x-ray imaging system 1 0 which can 

implement and perform the C-DRM processing algorithm of the invention. System 10 
includes a table 1 2 on which an object 1 4 is placed prior to being imaged (e.g., a 
human body). Object 14 is irradiated by x-rays emitted by an x-ray tube 1 6, driven by 
a generator 1 8 under exposure management control by a control and preprocessing 
computer 20. X-rays passing through the object 14 are attenuated in differing 
amounts depending upon the character of the body portions, and conditions therein, 
of the object 1 4 being imaged. Thus, the variously attenuated x-rays impinge upon 
and are detected by an image detector 22, the image signal output of which is 
O supplied to the control and preprocessing computer 20, all in conventional fashion. 

m The image detector 22 may be a digital detector or an analog detector. In use, an 

analog detector is coupled with a pick up tube or CCD resulting in a digital imager. As 
FU is likewise conventional, the preprocessing computer 20 controls various functions of 

r5 the x-ray imaging process including, e.g., offset, gain, bad pixel correction, exposure 

management (e.g., as to the x-ray tube 1 6 and generator 1 8) and scaling operations 
m and, further, preprocesses the received x-ray image data so as to produce a corrected 

™ x-ray image data output X to a display optimizing computer 30. 

LU 
□ 

U [001 6] A display-optimizing computer 30 includes a memory 32 for processing the 
image data X output by the preprocessing computer 20. Image data X is to be 
processed in accordance with the C-DRM algorithm of the invention by the computer 
30 and, further, performs related display optimizing processing, such as spatial and 
temporal filtering, zoom, and display control functions, for producing an optimized 
output signal Y. Signal Y is supplied to, and displayed by, display unit 34. Computer 
30 is programmed to perform functions described above and all exemplary 
embodiments described below, accordingly, as used herein, the term computer is not 
limited to just those integrated circuits referred to in the art as computers, but 
broadly refers to computers, processors, microcontrollers, microcomputers, 
programmable logic controllers, application specific integrated circuits, and other 
programmable circuits. 

[001 7] Figure 2 is a flow diagram of an EDR algorithm processor 40 of the prior art. An 
EDR function is computed and implemented in the following manner: 
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[0018] 



Y(i j) = Gamma[x(i j)-Boost[x(ij)]J (l) 
where: 

y(i, j) is the (i, j) th pixel value of the output image; 
x(i, j) is the (ij)* pixel value of the input image; and 

* (i, j), i.e., (x bar(ij)), is the local spatial mean intensity value of the 

(i j) * pixel, derived from a moving-average algorithm such as a BOXCAR average 
42. In practice, the x-rays are converted to photons which, in turn, are converted to 
electrons in an x-ray imager. The electrons are then digitized by an analog-to-digital 
converter ("A/D" or "ADC") and basic image corrections are performed, such as 
adjustments to gain, offset and scaling, after which the image is ready to be processed 
for display. 

The intensity value x of an input unit, such as a pixel, (i j) of an image 
is first processed by a BOXCAR function 42 that determines the local mean intensity 
value of the desired unit at that (ij) pixel location. (An V is used herein to designate 
an input intensity at a pixel location and, thus constitutes an individual scalar value; 
by contrast, an M X" designates the intensity image value, and thus is a vector value.) 
The BOXCAR function 42 utilizes a plurality of units, such as a neighborhood of 
pixels, which includes and is centered on the input pixel, to calculate the local spatial 
mean intensity value x (i j) ~ (see, terms of Equation (1), supra). 

As illustrated in Fig. 2, BOOST LUT 44 includes a look-up table 
(LUT) which specifies the intensity reduction of the input image signal x(i j) as a 
function of the local spatial mean intensity value J(ij). An adder (ADD) 46 
combines the (negative) output of BOOST LUT 44 and the input image signal x(i j) 
(see Equation (1)) and supplies the result to GAMMA LUT 48, which then 
compresses the result of the unsharp masking, or subtraction, operation of an ADD 46 
to 256 levels (8 bits per pixel, or 8 bpp). The LUTs 44 and 48 are indexed by the 
appropriate pixel intensity values given in equation (1). Thus, each of the BOOST 
LUT 44 and the GAMMA LUT 48 jointly manage both mean and contrast 
modification functions. 



[0019] 



Figure 3 is a flow diagram of a method 50 of implementing and performing the C- 
DRM algorithm in accordance with an exemplary embodiment of the invention. In one 
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aspect of the invention a method 50 is provided for compressing an intensity dynamic 
range of an input image to a reduced intensity dynamic range of an image display 
device. The method includes defining 52 a plurality of units of each frame of the input 
image and determining 54 a local mean estimate of the intensity of each of the 
plurality of units of the input image. The method further includes generating 56 a 
contrast modification processing value for each local mean estimate, and generating 
58 an output intensity value, for each local mean estimate, by using only the 
generated contrast modification processing value and an input intensity value. 

[0020] Fig. 4 illustrates a block diagram of a C-DRM processor 60 of the display- 
optimizing computer 1 2 of Fig. 1 , in accordance with an exemplary embodiment of 
r«5 the invention. Fig. 4 is labeled with bit size designations relevant to practical 

implementations of the apparatus implementing and performing method 50 of Fig 3. 
S3 Thus, image input X is designated "unsigned 1 2 bpp", i.e., an unsigned 1 2 bit per 

L}f pixel input image signal. Similarly, bit depths of intermediate processing results are 

Lj likewise labeled in bits per pixel ("bpp"). Further image input X may be encoded to 

* ' account for the framework of the specific invention. Encoding methods include, but 

y are not limited to, linear encoded data and square-root encoded data. 

gi [0021] Each block of the C-DRM processor 60 represents one or more functions, where: 

M [0022] 
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X 




input intensity value 


M x 




local mean estimate 


Cx 




contrast estimate 


FmO 




mean modification function 


FcO 




contrast modification function 


Y 




C-DRM output intensity value 


M Y 




modified local mean 


C y 




modified contrast 



The C-DRM processor 60 computes Y as follows: 



Cy = C x *f c {M x ) (3) 
M y =f u (M x ) (4) 



C ¥ Y 



Y = C r *M r =[C x *f c {M x ))*f M {M x ) (6) 



M 



LUT „ f.(M x )f m (M x ) 

M r 



(7) 



(8) 



Page 7 of 28 



# 



CombinerlnputValue = X* M M x)f m {M x ) 

My 



(9) 



Ul 
O 

w 



In accordance with an implementation of the exemplary embodiment of 
the invention shown in Fig. 4, the input intensity value X of an input pixel is linearly 
encoded. The image is first processed by a "BOXCAR" (or moving average) function 
62 to generate a local mean estimate M x . BOXCAR function 62 utilizes a plurality of 
units such as pixels, typically next-adjacent to and including the input pixel being 
processed, to calculate the local mean estimate M x . In an alternative embodiment, 
BOXCAR 62 may also use any reasonable array size to perform the BOXCAR 
average. BOXCAR function 62 acts as a low pass filter which attenuates or 
eliminates high-frequency components while leaving low frequency components 
untouched. A low pass filter separates structures of interest from variations in the 
detection output which are less relevant, e.g., gradually varying background. In 
experimental implementations, a 32x32 BOXCAR resulted in ringing along contrast 
filled vessels, heart to lung transition, heart to diaphragm, contrast filled catheter, etc. 
A 64x64 BOXCAR avoids that deficiency by spanning a sufficient spatial area, or 
extent, such that contrast-liquid filled coronaries do not unduly influence the mean 
estimate of the background and sudden background changes, such as heart to lung, are 
more gradual. BOXCAR filters of a particular size require specification of the image 
resolution as projected onto the panel to be defined. The scope of this invention is not 
limited to any particular BOXCAR size. 

The local mean estimate M x is then inputted to LUT1 64, which 
returns a contrast modification processing value (LUT IO ut) for each M x input in 
accordance with equations (2-8) above. More specifically, a contrast modification 
function f c (M x ) is multiplied by a mean modification function f M (M x ) to generate a 
first value. The first value is divided by the local mean estimate M x to generate a 
second value. A look-up table is accessed, with reference to the second value, to 
generate a contrast and mean modification processing value (LUTiotrr), also referred 
to herein as a contrast modification processing value. The contrast modification 
processing value is combined with the input intensity value X only in a multiplier 66. 



Page 8 of 28 





Multiplier 66 performs the multiplication function equation as denoted in equation (9) 
and outputs the combiner input value to a combiner/ex tractor 68. The bit 
extractor/combiner 68 extracts the integer portion of the output result of multiplier 66 
and discards the fractional component. In an exemplary embodiment, 
combiner/extractor 68 outputs a composite (unsigned) 12 bpp extractor output value 
to generate an output intensity value Y within the reduced intensity range. 



pixel is square-root encoded. The image is first processed by a "BOXCAR" (or 
moving average) function 62 to generate a local mean estimate M x . The local mean 
estimate Mx is then inputted to a LUT1 64, which returns a contrast modification 
processing value (LUT ]0 ut) for each Mx input in accordance with equation (10) 
above. More specifically, a contrast modification function fc(Mx) is multiplied by a 
mean modification function fivi(Mx) to generate a first value. The first value is 
divided by the local mean estimate M x to generate a second value. A square-root of 
the second value is determined, and LUT1 64 is accessed to generate a contrast 
modification processing value (LUTiout) which is combined with the square-root 
encoded input intensity value X in a multiplier 66. Multiplier 66 performs the 
multiplication function equation as denoted in equation (11) and outputs the combiner 
input value to a combiner/extractor 68. The bit extractor/combiner 68 extracts the 
integer portion of the output result of multiplier 66 and discards the fractional 
component such that the number of bits of the extractor output value correspond to a 
predetermined number of bits required by the image display device. In an exemplary 
embodiment, combiner/extractor 68 outputs composite (unsigned) 12 bpp extractor 




(10) 




(11) 




In an alternative embodiment, the input intensity value X of an input 
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output value to a LUT3 70. LUT3 70 performs a squaring operation on the extractor 
output value to generate an output intensity value Y within the reduced intensity range 
as denoted in equation (12). In an alternative embodiment, the extractor output value 
is not limited to 12 bpp, but rather any value that generates an output intensity value Y 
within the reduced intensity dynamic range such as, but not limited to, 16 bpp. 



[0023] 



Fig. 5 illustrates a block diagram of a C-DRM processor 80 of the display- 
optimizing computer 1 2 of Fig. 1 , in accordance with an exemplary embodiment of 
the invention. Fig. 5 is labeled with bit size designations relevant to practical 
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implementations of the apparatus implementing and performing method 50 of Fig 3. 
Thus, image input X is designated "unsigned 12 bpp", i.e., an unsigned 12 bit per 
pixel input image signal. Similarly, bit depths of intermediate processing results are 
likewise labeled in bits per pixel ("bpp"). 

[0024] Each block of the C-DRM processor 80 represents one or more functions, where: 

[0025] 
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X 

M x 

Cx 

F M () 

FcO 

Y 

M Y 
C Y 



input intensity value 

local mean estimate 

contrast estimate 
mean modification function 
contrast modification function 
C-DRM output intensity 
modified local mean 
modified contrast 



The C-DRM processor 80 computes Y as follows: 
M x 

Cy=C x *f c {M x ) 

C r '= Sigmoid (C r ) 
M r =f u {M x ) 

Y -My 



c, =- 



[0026] 




(13) 

(14) 

(15) 
(16) 

(17) 

(18) 



In accordance with the implementation of the exemplary embodiment of the 
invention shown in Fig. 5, an image is processed in accordance with equations (13- 
1 8) above. The input image is first processed by a morphological filter 82. 
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Morphological filter 82 is used to estimate the background from a current image. 
Morphological filter 82 is designed to determine an optimal mask from a given image 
in order to perform a digital subtracted angiography. In this embodiment, 
morphological filter 82 is used to estimate the image background for optimal 
compression with contrast preservation and generate an input image to the C-DRM 
filter. 

[0027] 
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More specifically, an input intensity value X of a pixel in the input 
image is first processed by the morphological filter 82 to generate a local mean 
estimate M x The local mean estimate M x is then inputted to an arithmetic unit 86, 
which processes X in accordance with equation (13) above, and outputs the contrast 
estimate C x to a multiplier 88. A contrast LUT 90 processes M* in accordance with 
the contrast modification function F c () of equation (14) above and supplies the 
contrast gain, the output f c (M x ) thereof, to multiplier 88 which performs the 
multiplication function of equation (14) and produces the modified contrast Cy output 
to a Sigmoid LUT 92 in accordance with equation (15) above. The sigmoid LUT 92 
function is applied to the modified contrast Cy output and used to enable linear 
preservation of small inputs and non-linear and asymmetrical compression of large 
inputs. Alternatively, a mathematical function may be used to replace Sigmoid LUT 
92. The use of the morphological filter for mean estimation and the sigmoid LUT for 
tuning contrast may be applied separately or jointly and are not dependent on each 
other. In another embodiment, Sigmoid LUT 92 can be used with a BOXCAR 
function, a morphological filter, or any other mean estimator such as, but not limited 
to, a median filter. 

Mean LUT function 94 processes the local mean estimate M x by the 
mean modification function F M () of equation (16), above, to output the modified local 
mean My. The mean modification function F M 0 compresses the low frequency 
component of the input intensity X, directly and separately from (i.e., independently 
of) the contrast component. A component 96 multiplies the Sigmoid LUT 92 output 
and Mean LUT 94 output to generate a value which is then combined with Mean LUT 
92 in an adder 98, outputting the C-DRM output intensity Y to display 34, thereby to 
reproduce the image thereon in the reduced intensity dynamic range supported by 
display 34 as depicted in Figure 1 . 



[0028] While the Invention has been described in terms of various specific embodiments, 
those skilled in the art will recognize that the invention can be practiced with 
modification within the spirit and scope of the claims. 
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